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Abstract. In this work we discuss a theory for entanglement generation, 
characterization and detection in fermionic two-particle interferometers at finite 
temperature. The motivation for our work is provided by the recent experiment by 
the Heiblum group, Neder et al, Nature 448, 333 (2007), realizing the two particle 
interferometer proposed by Samuelsson, Sukhorukov, and Biittiker, Phys. Rev. Lett. 
92, 026805 (2004). The experiment displayed a clear two-particle Aharonov-Bohm 
effect, however with an amplitude suppressed due to finite temperature and dephasing. 
This raised qualitative as well quantitative questions about entanglement production 
and detection in mesoscopic conductors at finite temperature. As a response to these 
questions, in our recent work, Samuelsson, Neder, and Biittiker, Phys. Rev. Lett. 
102, 106804 (2009) we presented a general theory for finite temperature entanglement 
in mesoscopic conductors. Applied to the two-particle interferometer we showed that 
the emitted two-particle state in the experiment was clearly entangled. Moreover, we 
demonstrated that the entanglement of the reduced two-particle state, reconstructed 
from measurements of average currents and current cross correlations, constitutes a 
lower bound to the entanglement of the emitted state. The present work provides an 
extended and more detailed discussion of these findings. 
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1. Introduction 

There is presently a strong interest in computation and information processing based 
on fundamental principles of quantum mechanics [Ij. Quantum information technology 
has the potential both to address problems that can not be solved by standard, classical 
information technology as well as to radically improve the performance of existing 
classical schemes. The prospect of scalability and integrability with conventional 
electronics makes solid state systems a likely future arena for quantum information 
processing. Of particular interest is the entanglement between the elementary charge 
carriers, quasiparticles, in meso- or nanoscopic solid state conductors. Entanglement, or 
quantum mechanical correlations, constitutes a resource for any quantum information 
process. Moreover, due to controllable system properties and coherent transport 
conditions, conductors on the meso and nano scale constitute ideal systems for the 
generation and detection of quasiparticle entanglement. This opens up for quantum 
bits based on the spin or orbital quantum states of individual electrons, the ultimate 
building blocks for solid state quantum information processing. 

To date quasiparticle entanglement has however remained experimentally elusive. 
In particular, there is no unambiguous experimental demonstration of entanglement 
between two spatially separated quasiparticles. A class of mesoscopic systems that 
appear promising for a successful entanglement experiment are conductors without 
direct interactions between the quasiparticles. It was shown by Beenakker et al [2] 
that fermions emitted from a thermal source can, in contrast to bosons, be entangled 
by scattering at a beam-splitter. This was originally discussed for electron-hole pairs 
[2] and shortly afterward for pairs of electrons [3], H] . Since then there has been a large 
number of works on entanglement of non-interacting particles, see e.g. [SI EJ [71 El [HI [ID] 
for a number of representative papers and also [11] for a review. 

Several of the entanglement proposals have been based on electrical analogs of 
optical interferometers and beam-splitter geometries. Such electronic systems are 
conveniently implemented in conductors in the quantum Hall regime, where electrons 
propagate along chiral edge states [121 US] and quantum point contacts constitute 
refiectionless beam-splitters [HI [151 [IB] with controllable transparency, see e.g. [T7] . 
Recent experimental progress on electronic Mach-Zehnder [181 [13 [23 [211 [22] and 
Hanbury Brown Twiss [23] interferometers has provided further motivation for a 
theoretical investigation of entanglement in such systems. In addition, the experimental 
realization [2l] of time-controlled single-electron emitters [251 126] in quantum Hall 
systems has opened up the possibility for a dynamical generation of entangled 
quasiparticles, entanglement on demand [271 [281 [29l [30]. 

In this work we will focus on the electronic two-particle, or Hanbury Brown 
Twiss, interferometer. A theoretical proposal for an implementation of this two-particle 
interferometer (2PI) in a conductor in the quantum Hall regime was proposed by two 
of us, P.S and M.B., together with E. V. Sukhorukov in Ref. [3]. Recently, the Heiblum 
group, including one of us, I.N., was able to realize the 2PI in a versatile system which 
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could be electrically tuned between with two independent Mach-Zehnder interferometers 
and a 2PI. In perfect agreement with the theoretical predictions [3], the two-particle 
interference pattern was visible in the current correlations but not in the average current. 
As discussed in Ref. [3], there is an intimate relation between two-particle interference 
and entanglement in the fermionic 2PI. Under ideal conditions, i.e. zero temperature and 
perfect coherence, two-particle interference implies that the two particle wave function 
is on the form 

l^.) = ^[|l)A|2)B-|2)A|l)i.]. (1) 

Here 1,2 denote the sources and A,B the sites of detection, as shown in Fig. [H The 
wavefunction l^^g) is maximally entangled, it is a singlet in the orbital, or pseudo spin, 
space {|1),|2)}. 

However, in the experiment ^J, ~ 25% visibility of the current correlation 
oscillations was observed. This indicates that both decoherence and finite temperature 
is important. Dephasing can qualitatively be accounted for [3ll|32l[33] by a suppression 
of the off-diagonal components of the density matrix |\E's)(\E's|- It was shown that at 
zero temperature the entanglement survives for arbitrary strong dephasing. The effect 
of finite temperature was not investigated at the time of the experiment. 

The experimental findings thus raised two important questions: are the electrons 
reaching the detectors at A and B entangled and if so, can this two-particle entanglement 
be unambiguously detected by measurements of currents and current correlators, the 
standard quantities accessible in electronic transport measurements? In our recent work 
[31] we provided a positive answer to both these questions. We first calculated the 
entanglement of the emitted two-particle state and found that the state was clearly 
entangled. Thereafter we showed that under very general conditions the entanglement 
of the reduced two-particle density matrix provides a lower bound for the entanglement 
of the emitted two-particle state. Since the reduced density matrix is possible to 
reconstruct tomographically by current and current correlation measurements [35] , this 
provides an unambiguous way to detect the entanglement of the emitted state. In the 
present paper we discuss these findings in more detail. 

2. The two-particle interferometer in optics and electronics 

Interference is most often investigated in structures that lead to a superposition of 
amplitudes of a single particle. However, in 1956, Hanbury Brown and Twiss (HBT) 
invented an optical interferometer based on correlations of light intensities [36l [37], an 
optical 2PI, see fig. [H The intensity interferometer allowed HBT to determine the 
angular diameter of a number of visual stars, not possible with available single particle, 
or Michelson, interferometers. The HBT intensity interferometer displays two distinct 
but fundamentally interrelated features: 
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• First, there is a direct statistical effect since photons from a thermal light source 
tend to bunch, whereas fermions would anti-bunch. This effect has been used in a large 
number of experiments in different fields of physics such as elementary particles |38j . 
solid state [HI [151 [16] cind free [39] electrons and recently cold atoms [40J. 

• Second, light from two different, completely uncorrelated sources gives rise to an 
interference effect in intensity correlations but not in the intensities themselves. This is 
the two-particle interference effect. In optics, various aspects of two-particle interference 
have been investigated extensively since the HBT-experiment, see e.g. [H] for a short 
review, and is still a subject of interest [12]. In electronics, only very recently was a 
fermionic two-particle interferometer realized [23], the subject of this work. 

Fundamentally both of these effects are related to the symmetry of the multiparticle 
wave function under exchange of two particles. We note that albeit the HBT-experiment 
could be explained by a classical electro-magnetic theory, a compelling quantum 
mechanical picture based on individual photons was put forth soon after the experiment 
[13]. Importantly, for fermions no classical theory exists. 




Figure 1. a) Schematic of the Hanbury Brown Twiss intensity interferometer used 
to measure the angular diameter of stars. Two uncorrelated points 1,2 on the star 
act as sources. The signal is detected at A and B. b) Schematic of the topologically 
equivalent two-particle interferometer (2PI) [3] with beam splitters C,D and biased, 
active (grounded, inactive) source contacts 1,2 (3,4). Detector regions A and B (red 
shaded) contain beam splitters and grounded contacts ±. 

To obtain a qualitative understanding of the physics of two-particle interferometers 
it is rewarding to compare the properties of optical, bosonic interferometers and 
electronic, fermionic interferometers. In Fig. [1] a schematic of a two-particle 
interferometer, topologically equivalent to the HBT-interferometer, is shown. A natural 
measure of the correlations between the particles at A and B is the probability to jointly 
detect one particle at A and one at B. An expression for this joint detection probability 
for photons was derived by Glauber [H]. In Ref. [3] this was adapted to detection of 
electrons. Here we consider the probability to detect one photon/electron in detector 
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Aa, a = ±, at time t and one in detector Bf3, /3 = ± at a time t + r, given by 

PAaB^ir) CC {bi^{t)b\^{t + T)bAa{t + r)bBp{t)) (2) 

The photon/electron creation operators at A are b]^^(t) = J dE exp{iEt / h)bl^^{E) , with 
b\^{E) creating a particle in Aa at energy E and similarly at B. For photons we consider 
thermal sources in 1 and 2 while 3 and 4 are left empty. A detector frequency window 
of size Auj is assumed, over which the distribution functions of the sources are constant, 
i.e. Auj ^ kT. For electrons we assume zero temperature and the sources 1 and 2 
biased at eV while sources 3 and 4 are grounded. Only quasiparticle excitations, E >0 
are considered. 

The probabilities are normalized such that J2a,p=±PAaBi3 = 1- Following the 
scattering theory for intensity/current correlations for bosons/fermions emitted from 
thermal sources [l5l HH] , we get 

PAaB^ir) OC \sAal\^\sBpi\^ [1 ± g{r)] + \sAa2\'^\sBl32\'^ [1 ± g{r)] 

I |2| |2 I |2| |2 

+ \SAal\ \SBf32\ + \SAa2\ \SBf3l\ 

± fi'(T) [sAal''^*Bl32''^BpiSAa2 + -SAal-S 5/32-55/31-5^02] (3) 

where g^r) = sir? {j jnTc) / {j /titcY contains the time dependence, with tq = h/eV 
the coherence time for electrons and 2/7rAc<j for photons. Here SAa2 is the amplitude 
to scatter from source 2 to detector Aa etc. The upper/lower signs ± correspond to 
electrons/photons. 

Several interesting conclusions can be drawn directly from Eq. ([3]): 

1) For r ^ Tc, (yf(r) approaches zero and PAaB/s is just proportional to the product 
of the two mean currents/intensities. The fermionic versus bosonic statistics of the 
particle plays no role. 

2) For shorter times, r < tq, g{T) is finite and the statistics is important. Note 
that, as pointed out above, that the statistics of the particles enter in two different ways. 

i) The first two terms in Eq. ([3]) describe a direct bunching (+) or anti-bunching (-) 
effect for two particles emitted from the same reservoir within a time r < tc- This 
effect would still be present if one of the sources 1 or 2 is removed. 

ii) The last two terms describe the two-particle, or exchange P5| B6]. interference, where 
the ± sign explicitly follows from the interchange of the two detected particles. This 
two particle interference is only present when both sources are active. 

For semitransparent beam-splitters A, B, C and D and coincident detection r <^tc 
we have 



PaolBP 



1 + a(3 cos I 

1 + COS( 



fermions , , 

bosons 



where is a scattering phase. From this expression a very important difference between 
bosonic and fermionic thermal sources is apparent: the visibility 

pmax JDmin 

pmax _|_ pmin \ / 

^AaBf3 "T ^AaBf} 
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of the oscillations is 1 for fermions but only 1/2 for bosons. This is directly related to 
the fact that while the emitted fermionic two-particle state is maximally entangled, the 
bosonic state is unentangled [19] . 

3. Fermionic two particle interferometer: theory 

In Ref. |3] we proposed an implementation of an electronic 2PI in a conductor in the 
quantum Hall regime, with electrons propagating along single, spin polarized edge states 
(see Fig. [2]). Two electronic reservoirs 1, 2 biased at eV act as sources for electrons while 
the reservoirs 3, 4 as well as the detector reservoirs are grounded. All reservoirs are kept 
at the same temperature T. Moreover, we consider here only the linear regime in voltage 
where electron-electron interactions can be neglected. This regime is relevant for the 
experiment [23]. The QPC's at A, B, C and D act as beamsplitters with transparencies 
Ta,Tb,Tc and respectively. 

a) , |Pn b) 



1 -c; 




B 



Figure 2. a) Fermionic 2PI implemented in a conductor in the quantum Hall regime, 
from 3]. See text for details, b) Schematic of the source part of the 2PI, with the 
the orbital states \1)a, |2)a, and \2)b for particles emitted out from the source 
towards the detectors are shown. 

The scattering amphtude sa+i = V^A-Rce*'^'*^, where Rc = I — Tc and (pAC 
is the scattering phase picked up by the electron up when traveling from C to A. 
Similar relations hold for the other scattering amplitudes. Note that the total phase 
= (pAc ~ 'Pad + <Pbd — <Pbc up to a constant term, given by 27r$/$o where $ is 
the magnetic flux threading the 2PI and $o = h/e, the single particle flux quanta. 
Importantly, the Corbino geometry in Fig. [2] with unidirectional edge states and 
reflectionless beam-splitters is topologically equivalent to the 2PI shown in Fig. [H 

3.1. Two particle Aharonov-Bohm effect 

The standard tools for investigating transport properties in mesoscopic electronic 
systems are average electrical current and current correlation measurements [17]. A 
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scattering theory calculation [l8] gives the average current at contact Aa 

(6) 

and similar at B(3. Here fv = 1/(1 + e^^'^^'^/^^) and / = 1/(1 + e^/'^^) are the Fermi 
distributions of the biased, 1,2 and the grounded, 3,4 reservoirs respectively. The 
irreducible zero frequency correlator 

SAaBp = J dt{AlAaiO)AlBpit)) (7) 

between currents lAa{t) = ^Aa + ^-^Aq(^) and I spit) = Isp + '^iB/sit) [IS] becomes 

SAaBf^ = ^JdE {K^.SBfSl + S\^,SBP2\') [ME) - f {E)f (8) 

These expressions are valid for arbitrary temperature but for the rest of the discussion 
in this section we only consider the zero temperature case. In particular, for the simplest 
possible case, with all beam-splitters semitransparent and energy-independent scattering 
amplitudes, we have 

lAa = I BP = SAaBp = ^ [1 + COS (f)] (9) 

While the average current is a function of QPC-transparencies only, the current cross 
correlator depends also on the phase 0. Since this phase is proportional to the magnetic 
flux $ threading the 2PI, we call this a two-particle Aharonov-Bohm (AB) effect. 

Interestingly, we can directly relate the coincident detection probability in Eq. ([3]) 
at times r <^ tq with the currents in Eq. ([6]) and the zero frequency noise correlators 
in Eq. (JED as [^(0) = 1] 

PAaBpi^) OC \sAalS*Bpi + SAa2S*Bp2\^ + (^AalP + |'5yla2p) (k_B/31P + k_B/32p 
OC SAaBp + 2TclAalB/3 (10) 

This is a direct consequence of fermionic anti-bunching, leading to a filled stream of 
electrons emitted from the source reservoirs and hence making long time observables an 
effective average of many individual, short time, single and two-particle events. 



3.2. Entanglement 

The connection between this two-particle Aharonov-Bohm effect and entanglement can 
be seen by considering the many-body ground state {"^m) of the electrons injected into 
the 2PI. Electrons at different energies are independent and the many-body state at 
zero temperature is thus a product state in energy 

|^.„)= n a\{E)al{E)\0) (11) 

0<E<eV 

where |0) is the filled Fermi sea and a\{E) creates an electron at energy E, incident from 
reservoir 1. Adopting the formalism of Ref. |^ we first define |\['m(i?)) = a\{E)a2{E)\0) 
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the injected state at energy E. We have the scattering relations at the two source beam 
sphtters, suppressing energy notation 

/ . X / . ^^^^ 

for incoming (a's) and outgoing (b's) electrons. The primed scattering amplitudes thus 
describes particles incoming from the unbiased sources. This gives the emitted state 
for the electrons at energy E, after beam-splitters C, D but before impinging on the 
detector beam splitters A, B, as 

\^out{E)) = (rcb\, + tM (rDb\, + tz^&L) |0) (13) 



Since we are interested in entanglement between particles in the two, spatially separated 
detector regions A and B we project out the part of the wave function with one particle 
in A and one in B yielding the normalized wavefunction 

l^ABiE)) = {rctDb^A, - rntcb^Ai) |0) (14) 



with = |r£)tcP + kc^Dp = RcTd + RdTq the normalization constant. Here we 
introduced the transmission and reflection probabilities of the source beam splitters as 
Tc = |tcP = Wc? and Rc = |rcp = |r^P = 1 - Tc for C and similarly for D. To 
make this more transparent we can, since the two particles live in well separated Hilbert 
spaces, introduce the Dirac notation = &^i|0) etc, and write 

l^ABiE)) = 4tt Md|1)a|2)b - tcrc\2)A\l)B] (15) 



which for semi-transparent beam splitters (and scattering phase = 0) reduces to the 
singlet state |\l/s) in Eq. ([1]). The orbital states are shown in Fig. [2] 

The entanglement of the state |\I^^b(-E)) can conveniently be quantified in 
terms of the concurrence C [50], which ranges from zero for an unentangled state 
to unity for a maximally entangled state. Working in the computational basis 
{\1)a\1)b, \1)a\2)b, \2)a\1)b, \2)a\2)b}, for the pure state |^ab) in Eq. ([HD we have 

C=\{^AB\icry^ayWA^)\ (16) 

where 1^^^) is I^E'ab) with all coefficients complex conjugated, ay a Pauli matrix and 
the direct, tensor product. We thus find for |\1/ab) the concurrence 
2 , ,2 



C = -IrctcrDtol = RcTcRdTd (17) 

which reaches unity for semitransparent beam splitters, i.e. for the singlet state in Eq. 
([I]). Note that the normalization factor A^ is maximal, equal to 1/2, for semitransparent 
beam splitters. This demonstrates that at most only half of the particles injected from 
1 and 2 lead to split pairs, with one particle emitted towards A and one towards S, i.e. 
a maximal pair emission rate of 1/2. For a measurement during a time r the maximum 
concurrence production [H] is thus A/'/2, where N = reV/ h the number of pairs injected 
from 1 and 2 in the time r and energy interval < -E < eV 
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3.3. Dephasing 

There are several microscopic mechanisms that can lead to dephasing, typically 
suppressing the two-particle interference. For low temperatures it is commonly believed 
that the dominatinating mechanism for dephasing is electron-electron interactions, but 
this is still a topic of ongoing research and goes beyond the scope of the present work. 
Here we consider no specific mechanism but model dephasing qualitatively by coupling 
one of the interferometer arms to a dephasing voltage probe [5T1 [52| [53l IM] . In this 
context we point out a recent experiment [55j: a voltage probe was coupled, via a tunable 
quantum point contact, to one arm of a Mach Zehnder interferometer in the quantum 
Hall regime, demonstrating controllable dephasing. Considering semitransparent beam 
splitters, the dephasing probe coupled with a strength < 7 < 1 lead to a modification 
of the current correlator in Eq. ([9]) to [56] 

^t% = ^[l + 7«/5cos0] (18) 

From this expression it is clear that 7 enters as a decoherence parameter; decreasing 7 
from 1 to leads to a suppression the phase dependence of the current correlator. In the 
presence of dephasing the emitted state is no longer a pure state, it is instead a mixed 
state described by a density matrix aAB- Considering zero temperature, working in the 
computational basis the result for S'^^^^ corresponds to a suppression of the off-diagonal 
components of \'^ab){'^ab\ o'ab as 

/ \ 

1 -7 



1 

2 -7 1 
V / 

The concurrence for a mixed state is EO 



(19) 



C = max|^Ai - ^As - V^3 - V^4,0| (20) 

where Aj, z = 1 — 4, are the eigenvalues in decreasing order of o'AB{o'y ^ o'y)o'ABi^y ® ^y) ■ 
We then have 

C = i (21) 

This means that the entanglement persists even for very strong dephasing [HH [32l [33] . 
This is a consequence of the 2PI-geometry, where scattering between the arms, i.e. 
pseudo spin-flip scattering, is prohibited. 



3.4. Fermionic two particle interferometer: experiment 

Very recently the electronic 2PI was realized experimentally by Neder et al. In the 
experiment, in the quantum Hall regime, it was possible to electrically tune the system 
between two individual Mach Zehnder interferometers and a 2PI, as shown schematically 
in fig. [3l The authors first tuned the system to two Mach-Zehnder interferometers and 
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measured the single particle interference in the average current for each interferometer. 
They found a very large visibility in both interferometers, around 80%. They also 
determined the periods of the single particle AB-oscillations as a function of both the 
area and the magnetic flux enclosed by the interferometers. Thereafter the system was 
tuned to a single 2PI. As predicted by theory [3] no single-particle AB-oscillations in 
the average current were observed but the current cross correlations displayed clear 
two-particle AB-oscillations, with an amplitude 25% of the predicted coherent, zero 
temperature value. By measuring also the period of the two-particle oscillations as 
a function of interferometer area and enclosed flux and comparing to the sum of the 
periods for the two Mach Zehnder interferometers, the two-particle nature of the AB- 
oscillations could be established beyond doubt. 

In the experiment semitransparent beam splitters were used, Tc = = 1/2. For 
the current cross correlations, theory for finite temperature and dephasing [56] predicts, 
for A+,B+, 



a) 




b) 




Figure 3. Fermionic two-particle interferometer implemented in a conductor in the 
quantum Hall regime in Ref ^23J. a) Figure reproduced from Ref. |23j . Micrograph 
of the sample, b) Left: The system in the two Mach Zehnder interferometers 
configuration. Right: The system in the 2PI configuration. 




(22) 



The temperature dependence is fully contained in 




(23) 
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Figure 4. Figure reproduced from Ref. [23]. Experimental demonstration of the two- 
particle AB-effect. Current cross correlation displaying clear oscillations as a function 
of the effective interferometer area and enclosed magnetic flux. 

varying from unity for kT <^ eV to zero for kT ^ eV. The effect of finite temperature is 
thus to suppress the overall amplitude of the current cross correlation oscillations. In the 
experiment, the applied bias was 7.8^V. The electron temperature was estimated from 
independent auto-correlation measurements to be lOmK. This yields the temperature 
suppression factor H = 0.78. A direct comparison to Eq. (1221) then gives the oscillation 
amplitude H'j = 0.25, i.e. 7 = 0.32, a substantial dephasing. 

4. Finite temperature state 

Our main aim of this work is to theoretically investigate the effects of finite temperature 
on the entanglement of the state emitted out from the source, towards the detectors. 
A prerequisite is to obtain both a qualitative and a quantitative description of the 
emitted many-body state at finite temperature. We consider the experimentally relevant 
situation with all source and detector reservoirs kept at the same temperature T. Due 
to the finite temperature, not only the electrons emitted from the source in the energy 
range < E < eV are of interest, we must in principle take into account particles 
emitted from all reservoirs at all possible energies. However, due to the chiral geometry 
of the 2PI in Fig. [21 particles emitted from the detectors can never scatter back to the 
detectors, i.e. detector cross talk is topologically prohibited. The particles arriving at 
the detectors thus all originate from the source reservoirs and we can focus on the many 
body state emitted by source 1 to 4. We note that in the slightly different geometry 
realized experimentally [23], there is the possibility for scattering between the detectors. 
It can however be shown [5^ that this does not influence the entanglement of the emitted 
state. 

At finite temperature the state injected from the sources is mixed and described by 
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a density matrix [llj 

Pin = YlPiniE) 

E 

Pin{E)= n [[l-/.(i?)]|0)(0|+/4E)4(£;)|0)(0|a.(E)] (24) 
ft=i 

where fniE) is the Fermi distribution of source reservoir k = 1 — 4. The outgoing state 
is then obtained by inserting the scattering relations of Eq. f|T2l) int Eq. flM|) . 

One can see from Eq. flM|) that the effect of finite temperature is to give rise to 
states with to 4 particles emitted at a given energy. For the terms of interest, i.e. with 
at least one particle at both A and B, there is at finite temperature the possibility for 
e.g. two particles at A and one at B etc. These terms are of central importance in the 
discussion below. 



5. Projected two-particle density matrix 



A theory for entanglement production in non-interacting ^ conductors at finite 
temperature was presented by Beenakker [TT] and along similar lines in closed condensed 
matter systems by Dowling, Doherty and Wiseman [58]. At a given energy, only the 
component of the emitted many-body state with one particle in detector region A and 
one in B has nonzero entanglement. Moreover, as emphasized in Ref. [58], only this term 
describes two particles which each live in a well defined 2x2 Hilbert spaces at A and 
B respectively, i.e. two coupled orbital qubits. We point out that this definition does 
not take into account occupation- number, or Fock-space entanglement. The first step 
is thus to project out the two-particle component from the many-body wave function, 
which is accomplished by the projection operator 



n 



Ua^Ub, Ha =^qi(1 -"-02) +'^a2(l -^qi) (25) 

where uaj = b^jbAj with j = 1,2 etc is the number operator (suppressing energy 
notation). This yields the projected density matrix 

Pp{E) = Up{E)U (26) 

The elements of the density matrix Pp{E) are conveniently calculated from the relation 
[581 



[pp{E)]ij,ki = {Ub^ibijbBkbAi'n) (27) 

where, for any operator X, {X) = tr[Xp] is the standard quantum-statistical average. 
Some algebra gives the projected density matrix, formally equivalent to the density 
matrix calculated in [TT], Eqs. (B9) - (B13), 

fx 



PpiE) 



fffl 





V 



.12 

-21 



.21 
-21 







\ 


X J 



i2t 
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where x = ^ eV/kT ^^^^^ j ^^^^ ^^le Fermi distribution functions of the grounded and 
biased source reservoirs respectively. The coefficients 

c\1 = {Rc[1-x] + x){Td[1-x] + x), 
cll = {Tc[l-x]+x)iRD[l-x] + x), 



cll = {cllr = -l^RcTcRDTDe'^''{l - X? (29) 
with 00 an overall scattering phase of the beam splitters C and D. Thus, only the prefac- 
tor — fY depends on energy. As for the zero temperature case we have introduced 
dephasing as a suppression of the off-diagonal components of the density matrix. It 
follows from Eq. fl2Sl) that finite temperature leads to 

i) an overall modification of the energy-dependent probability for two-particle emission 
via the pref actor (1 — fY fy. 

ii) a suppression ~ {1 — xY of the off-diagonal components, equivalent to the effect of 
dephasing. 

iii) a finite amplitude for the diagonal density matrix elements [pp(£')]ii^ii and 
[pp(£')]22,22, i-e for two particles being emitted from either sources 1,3 or 2,4. 

Additional insight follows from writing the projected density matrix as 

p,{E) = (1 - ffpy [xpf''' + (1 - (30) 
where the diagonal density matrix 

pMag = ^ i + (1 _ ^)[^^ ^ i + i ^ (31) 

with the zero temperature single particle density matrices pA = -Rc|l)(l| +-R_d|2)(2| and 
Pb = Tc|l)(l| +Td|2)(2|. The density matrix 

p^"* = RcTd\12){21\ + RdTc\21){12\ 



- 7V7ci?c7^Di?D[e**°|21)(21| + e-^<^«|12)(12|] (32) 

results from the two-particle interference. Here we used the shorthand notation 
|12) = |l)yi|2)B with (21| = (112))"'' etc. Note that the effect of decoherence enters 
as a suppression of the two-particle interference [v]/*"*) — p*"*^ where = 

^/R^\l'2) - e''^'>VTcRE\21). 

Writing Pp{E) in the form in Eq. fl30l) shows that, taken the energy dependent 
prefactor fy{l — f'Y aside, the effects of finite temperature can be viewed as follows: 
First, the amplitude of the two-particle interference component p*"* is suppressed with 
increasing temperature as ~ {1 ~ xY- Second, the density matrix acquires a purely 
diagonal component Pp*"^ with an amplitude ~ x (note that tr[pp*"^] = 4, independent 
on temperature). 

For the entanglement, following [llj we introduce o"p and Wp{E), the normalized 
density matrix and the emission probability of the emitted two-particle state 
respectively, defined from 

Pp{E) = Wp{E)(Tp, 

Wp{E) = tT[pp{E)] = (1 - fYfy[{RcTD + TcRdYI - xY + 4x] (33) 
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where we note that CTp is independent on energy. The emission probabihty Wp{E) is thus 
the probabihty, per unit energy, that the (normahzed) two-particle state is emitted. 
The concurrence production per unit energy is then 



X ma.|4,VflcTcfloT„ - -p(j^,0| (34) 

and the total entanglement production during a time r, Cp = (r/h) J dECp{E), is then 

(^f = Tcv/h) 



C, = —ma. |4Vrcflcrofl. - -pp^. 0| . (35) 

We denote this the projected entanglement. As shown in Fig. [5l Cp decreases 
monotonically as a function of T. It reaches zero at a critical temperature given 
by 




.T^ = e^l„|xi±i2^^^1 (36) 

A-i^RcTcRnTD - 




Figure 5. a) Entanglement production Cp/M (blue, transparent) and Cr/M 
(green, opaque) as functions of temperature kT/eV and coherence 7 for the semi- 
transparent 2PI. b) Parameter Q as a function of kT/eV (blue line). Values 
0.25, 10~^, 10"^, 10"'^, 10""' shown (gray lines). Figure reproduced from Ref. 

For semi-transparent beam-splitters and zero dephasing, 7 = 1, the entanglement thus 
survives up to [TT] kT^ = 0.57eV. 

Inserting the parameter values from the experiment, we get Cp ~ 0.1 A/" and 
C((Tp) ~ 0.3, i.e. the state emitted by the 2PI is clearly entangled. Importantly, the 
effect of finite temperature is essentially negligible, the reduction in entanglement comes 
from decoherence. 

The entanglement of the projected density matrix is the entanglement one could 
access, had one been able to do arbitrary local operations and classical communication 
between A and B, i.e. fully energy and particle resolved measurements. Under realistic 
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conditions this is not possible, the accessible physical quantities are currents and current 
cross correlators. Is it possible to determine the projected entanglement with such 
measurements? The answer to this question is no, for two main reasons: 

i) As discussed above, at nonzero temperatures it is not only the biased source reservoirs 
which emit particles but also the grounded source reservoirs do. As a consequence, 
there is a finite amplitude for emitted states with two-particles at A and/or at B. These 
unentangled states contribute to currents and current correlators, which results in a 
detectable state with suppressed entanglement. 

ii) The current and current correlators provide information on the energy integrated 
properties of the many-body state, not on the emitted state at each energy. This 
lack of energy-resolved information leads to a further suppression of the detectable 
entanglement. 

Clearly, these effects of the thermally excited Fermi sea constitute generic problems 
when trying to detect entanglement in mesoscopic conductors. 

As a remedy for these finite temperature read-out problems it was suggested 
to work with detectors at very low temperatures [llj. Another idea was recently 
presented by Hannes and Titov [8] . They investigated detection of entanglement at finite 
temperatures via a Bell inequality and proposed to introduce energy filters at the drains. 
However, both schemes [HI [8] would lead to additional experimental complications 
in systems which already are experimentally very challenging. Our idea is instead to 
investigate what information about the projected entanglement can actually be deduced 
from current and current correlation measurements. 

In this context we also mention the recent proposal by Kindermann to produce 
and detect entangled electron-hole pairs in graphene via a Bell inequality formulated 
in terms of the transport part of the current cross correlators [l6], i.e. by subtracting 
away the thermal equilibrium correlators from the finite bias ones. In our work [34] we 
proposed a similar scheme for a general mesoscopic conductor. However, as was pointed 
out in [21] and is further discussed below, it is important that one performs a detailed 
comparison of the projected entanglement and the entanglement obtained from current 
cross correlation measurements. Without such a comparison, there is the possibility that 
one concludes, based on correlation measurements, finite entanglement where there is 
none, i.e. the projected entanglement is zero. 

6. Reduced two-particle density matrix 

We first consider the expression for the current and zero frequency current cross 
correlators at contacts A-\- and B-\- at finite temperatures. We have 




(37) 
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where (An^+AnB_|_) = {nA+nB+) — {nA+){nB+) is the irreducible correlator. As 
discussed above, the many-body state incident on the detectors originates from the 
sources. It is the properties of this state that determines the observables {nA+), {^13+) 
and {AnA+AnB+) and thus establishes a connection between the emitted state and the 
physical quantities accessible in a measurement. 



6.1. Energy resolved reduced density matrix 

In order to better understand the readout problem discussed above, we first discuss the 
energy resolved properties of the emitted state. If one would have access to energy filters, 
as proposed in [8], or would be working at zero temperature, by combining current and 
current cross correlations it would be possible to get direct access to the energy resolved 
quantities {nA+), {nB+) and {AnA+AnB+). As is discussed below, by a suitable set of 
measurements with different settings of the beam splitters at A and B one could then 
tomographically reconstruct the (unnormalized) density matrix of the state emitted out 
from the source beam splitters C and D, pf, with elements given by 

[pfh,M = {b\fii^bBkbAi) (38) 

We denote the energy resolved reduced density matrix. 

By comparing pf with the expression for the projected density matrix in Eq. fl28l) 
we see that it differs by the projection operators. Consequently, the reduced density 
matrix contains also the contributions from processes with more than one particle at A 
and/or at B. After some algebra we find the density matrix 



Pf = (1 - fffl 



/ X \ 

c}2 dl 

4? ci\ 



(39) 



V X J 

where we introduced x = x/[(l ~ /v)(l ~ /)] and the coefficients 



^12 



{Rc[l-x] + x){Tn[l-x] + x), 

{Tc[1-x]+x){Rd[1-x] + x)- (40) 

A comparison to the projected density matrix in Eq. ( 128|) shows that pf only differs 
formally from Pp{E) by the change x X at a number of places. This has the 
consequence that the normalized density matrix o".f = pf /w^, with = tr[pf] depend 
on energy. That is, in contrast to pp both the normalized, emitted two-particle state 
as well as the emission probability depend on energy. Quahtatively, as discussed above, 
the difference between p^ and Pp{E) arises from the fact that also states with more than 
one particle at A and/or B contribute to pf but not to Pp{E). Writing pf on a form 
similar to Eq. ( |30l) one sees that these three and four particle states contribute only to 
the diagonal part of pf . 
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Turning to the entanglement, the concurrence production = w^C{af) at energy 
E is then 

rE (1 - x?B{l - f? 
- 2 

X max [A.^RcTcRnT. - ^^^.^^^^^^^ _ _ ^y ^] (41) 

From the expression for the concurrence it becomes clear that the separable three and 
four-particle states are detrimental for the entanglement. Hence, finite temperature 
leads to a stronger suppression of the reduced, energy resolved density matrix than of 
the projected one. This is illustrated in fig. [H] where the corresponding concurrencies 
are plotted for semitransparent beam-splitters and different values of kT/eV. As is 




Figure 6. A comparison of tlie concurrence production rates (dashed) and Cp{E) 
(solid), as a function of energy for Tc = = 1/2 and different ratios eV/kT. 



clear from the figure, there is an energy Eq above which the concurrence is finite (up to 
E — >• oo). The energy £"0 is given by the condition C^{Eq) = 0, as 



Eq 



kT (^\n[2]-\n (1 - x)^! + A^RcTcRdTd - (1 + x) 



(42) 



What is moreover clear from Fig. [6] is that, for all energies, Cf{E) < Cp{E). The 
difference is obvious for energies E < Eq, where = 0. At these energies the 
probability for emission of separable three and four particle states is thus large enough 
to completely suppress the entanglement of the reduced density matrix. 

Importantly, the relation Cf{E) < Cp{E) holds for all settings of the beam splitters 
Tc and To, as is clear by comparing Eqs. ( IMll and ( HTl) . The reason for this is that 
the reduced density matrix contains contributions from all individual particle density 
matrices 0"^ with z,j > 1 (e.g. (J12 describes one particle at A and two at B) while 
the projected density matrix only depends on an. Since all 0"i2, cr2i, (T22 are separable 
and the concurrence is a convex quantity, i.e. C{j>iai + P20'2) < PiC{ai) +P2C(c"2) 
for pi + P2 = 1, the concurrence is always smaller than Cp{E). We point out that 



Entanglement at finite temperatures in the electronic two-particle interferometer 18 



this carries over to the total concurrence production found by integrating Eq. fl41l) over 
energy (resuh not presented here). 

It follows from Eq. (H2i) that for a critical temperature TJ^ the energy Eq ^ oo, 
i.e. the entanglement is zero for any energy. Interestingly, this happens for the same 
temperature as for the projected concurrence, Eq. ( |36ll . 



6.2. Finite temperature reduced density matrix 

Importantly, at finite temperature, without any energy filters, we do not have access to 
the energy resolved quantities discussed above, only to the total currents and current 
correlators measured at contacts Aa, Bj3. In Ref. [35] it was discussed how to, at zero 
temperature, tomographically reconstruct the reduced density matrix using currents 
and current correlations. Extending this scheme to nonzero temperatures it is natural 
to define the finite temperature reduced density matrix via the relation 



{Ve'^/hf 2Veyh ~ 



lAa®lB(^]pr}- (43) 



RrTn rfi? 



d^l RnTn 



(44) 



We emphasize that pr is reconstructed from observables already integrated over energy 
and does hence not depend on energy. Also note that pr is not given by integrating 
over energy, in fact the difference between the two density matrices is further discussed 
below. 

In Eq. US] the orbital current operators in the local basis {|1), |2)}, including the 
rotations at the detector splitters, are I^^ = (1 + ohia ■ cr)/2 and I^p = (1 + (^t^b ■ o")/2, 
with Ha ■ d' = SaO'zSa and ■ cr = SbO'zS^b where a = [o'x,0'y,az] a vector of Pauli 
matrices and 5*^^ (5*^) the scattering matrix of the beam splitter at A (B). 

Making use of the results for finite temperature current and current correlations in 
[56] we obtain the reduced density matrix 

/ RcTcil -H) \ 

Pr = 

V RdTd{1-H)J 

where d^ = {dH)* = -H'^^/RJI\^Rjr^e''t'^ . Comparing pr to both pp{E) and pf in 
Eqs. (l28l) and (l39l) it is clear that the qualitative effect of finite temperature is the same 
for the reduced density matrix. The quantitative effects are however different. First, the 
temperature dependence enters via H rather than via Xi giving a much stronger effect 
of finite temperature. This is the effect of having access to energy integrated quantities 
only. Second, in the expression for the average current in Eq. (!37ll . in the integrand 
one subtracts / which arises due to particles flowing out of the detector reservoirs. This 
yields smaller diagonal terms, to be further discussed below. 

It is illuminating, just as for Pp{E), to write pr as a sum of a diagonal and an 
interference part, 

Pr = {l-H)[pA(»PB]+Hp''''K (45) 



IC-l-D "21 

ill RdTc 
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From this we see that the effect of increasing temperature is to monotonically increase 
the amphtude for the separable product state pa ® Pb, while the amplitude of the 
interference component is suppressed. We can thus conclude the following properties 
for all three density matrices Pp{E),pf and pr- 

i) At zero temperature they all reduce to the same expression, p*"*. 

ii) Increasing temperature leads to a monotonic suppression of the two-particle 
interference component. 

iii) Finite temperature introduces an additional diagonal component, different for the 
three density matrices. 

Turning to entanglement, introducing the normalized reduced density matrix ar we 
can write 

Pr = WrCr 

Wr = tT[pr] = [RcTc + RdTd] ' H) + RcTn + RdTc- (46) 

We then define the total entanglement production during a time r as Cr = MwrCi^ar). 
It is 



Cr = 2Afmax{^TcRcTDRD[Hil + 7) - 1], 0} (47) 

here called the reduced entanglement. As Cp, decreases monotonically with increasing 
T. It reaches zero at a critical temperature given by the relation 

HiTl) = (48) 
1 + 7 

For perfect coherence, 7 = 1, we have /cTJ = 0.28eV^, close to one half of kT^. 
Importantly, in contrast to T^, is independent on the setting of the beam splitters. 

By comparing the expressions for the two quantities of main interest, the projected 
and reduced concurrencies, Cp in Eq. (135|) and Cr in Eq. (147|) . we can conclude the 
following: 

i) For both Cp and Cr the origin of the entanglement is the two-particle interference, in 
fact the component p™* gives rise to the positive term 2^^H'y^/TcRcTDRD, identical 
for Cp and Cr- 

ii) For both Cp and Cr finite temperature introduces a negative term, 
-AfH/[2smh\eV/2kT)] for Cp and -2Af{l - H)^TcRcTdRd for C,, which leads 
to a suppression of the concurrence. These terms arise from the separable, diagonal 
components of the corresponding density matrices. 

7. Entanglement bound 

Comparing Eqs. ( l35l) and ( H71) quantitatively we find that Cp > Cr for 

QiT) = — -^-T-, -, 7 < JTcRcTdRd. (49) 

independent on 7 (see Fig. [5]). Consequently, for beam splitters away from the strongly 
asymmetrical (tunneling) limit, the reduced entanglement constitutes a lower hound for 
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the projected entanglement. In the tunneling limit, however, the reduced entanglement is 
larger than the projected one. Thus, in contrast to the energy-resolved reduced density 
matrix pf, pr can be more entangled than pp. The origin of this difference is, as pointed 
out above, that when calculating (and measuring) pr the average currents flowing out 
from the detector reservoirs are subtracted, yielding a smaller diagonal component and 
hence a larger entanglement Cr- Importantly, since the transparencies Tc and can 
be controlled and measured via average currents in the experiment, it is always possible 
to verify independently that the condition in Eq. ( H9|) is satisfied. 

Turning to the experiment [23], for the relevant parameters we have Q{T) ^ 
4 X 10""^ <ti \/ RcTcRdTd ~ 0.25, showing the validity of the bound. However, 
Cr ~ O.OIA/" and based on the measurement [23] no conclusive statement can be made 
about Cr and hence not about Cp. In order to detect entanglement via measurements 
of currents and current correlations, one thus need to work at even lower temperature 
and further reduce the dephasing in the experiment. 

A more detailed understanding of this finite temperature readout problem can be 
obtained by comparing the properties of (Xp and 0"^.. For perfect coherence 7 = 1 and 
identical beam splitters Tq = = T = 1 — IZ one can (up to a local phase rotation) 
write 

(Tp/r = ^^p/A ® i + (1 - ep/.)|^.)(*.| (50) 

a Werner state [59j, with singlet weight [I^E's) is the singlet in Eq. ([T])] 

^ 2nTsmh\2eV/kT) , _ , _ JL_ 

l + 27^Tsinh2(2eV/A;T' ^'^ 2-H ^ 

Increasing kT/eV from zero, ~ 2e~^^^/'^^/(7?.T) becomes exponentially small while 
^r ~ kT/eV increases linearly. These qualitatively different behaviors, clearly illustrated 
in Fig. [5], are a striking signature of how a small kT/eV, having negligible effect on 
C((Tp), leads to a large suppression of C{ar)- 

From Eqs. (l35l) and (H71) follows also a counter- intuitive result: finite amplitude of 
the AB-oscillations is no guarantee for finite two-particle entanglement. This is apparent 
for cTj, in the limit of no decoherence 7 = 1 and identical beam splitters Tc = T^, since 
a separable Werner state, ^r > 2/3, can be decomposed [50] as 

^. = ^El0n)(0nl®l</>n)(0fl (52) 
^ n=l 

with the normalized states at A and B 

10^)= cose;^/^|l) + e-[^-^"l/43i^^A/B|2^^ 

ef^ = et"" = atan[//^], 0^^^ = Of"" = -acot[//^] 

= ^^^/#^ ^ +(-)forA(5) (53) 

This classically correlated state gives, via Eq. fj43|) . AB-oscillations with amplitude 
2(1 — ^r)/(2 — ^r) = H. Moreover, the reduced local single particle states are completely 
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featureless, ti sicTr) 



i/2 which means that there is no single particle 



Aharonov-Bohm effect. The existence of classically correlated two-particle states giving 
rise to Aharonov-Bohm oscillations in the current cross correlations but not in the 
currents provides further motivation for a complete tomographic reconstruction of the 
reduced density matrix in order to provide an unambiguous experimental demonstration 
of entanglement. 

8. Detecting entanglement: Quantum State Tomography and Bell 
Inequality 

8.1. Quantum state tomography 

As pointed out at several places above, the reduced density matrix can be reconstructed 
by a suitable set of current and current correlations measurements with different settings 
of the beam sphtters parameters, i.e. different n^,nB. A detailed description of this 
scheme is given in [35]. Here we only emphasize that the necessary tools, controllable 
reflectionless electronic beams splitters and phase gates, are experimentally available, 
as demonstrated in e.g. [HI [191 [20l [2ll [221 [23] 

8.2. Bell Inequality 

Another widely discussed [611 [311 [21 [3l [5l [6] approach to detect the entanglement in 
mesoscopic conductors is to use a Bell inequahty. Violation of a CHSH-Bell inequality 
[62] formulated in terms of currents and low-frequency current correlations demonstrates 
finite entanglement of pr- We point out that an optimal Bell test, requiring control over 
all three components of ua and ub, demands the same number of measurement and 
level of experimental complexity as a tomographic reconstruction of pr- The CHSH-Bell 
inequahty is 



where flsp/r is the Bell parameter for the projected/reduced state. The Bell parameter is 
formally determined by the projected/reduced density matrix ap/r and different settings 
of the detector beam splitters, reaching its maximum value ^^pj^ for an optimal setting 
of ua and ub. From ap and (Jr above, we can, using Ref. [63], calculate the maximal 
Bell parameters. For symmetric beam splitters, Tc = Td = T we have the simple result 



where the singlet weights 1 — and 1 — are given in Eq. ( ISTI) . This shows that the 
effects of decoherence and finite temperature enters separately in the Bell parameter. 
Moreover, as pointed out in Refs. [311 [321 [33] , at zero temperature a Bell inequality 
can in principle be violated for arbitrary dephasing. We also point out that a detailed 
investigation of conditions for violation of a Bell inequality in the presence of dephasing, 
in the solid state, was recently performed in Ref. |65j . 




(54) 




(55) 
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The limiting value for violation ^^^p/j, = 2 for T = 1/2 plotted in Fig. El It is 
clear that for the values kT/eV and 7 of the 2PI-experiment, while flsp < 2 in principle 
can be violated, a detection of entanglement by violating i^^,. < 2 is not possible. This 
demonstrates in a striking way the known fact [591, [64j that there are entangled states 
that do not give a violation of a Bell Inequality. 

9. Conclusions 

In conclusion, we have investigate the effect of finite temperature on the entanglement 
production and detection in the fermionic two-particle interferometer, presenting an 
extended discussion of the results in Ref. [3l]. A calculation of the entanglement of the 
two-particle state projected out from the emitted, finite temperature many body state 
shows that the state emitted in the two-particle interferometer in the experiment by 
Neder et al [23] is clearly entangled. By comparing the entanglement of the projected 
two-particle state with the entanglement of the reduced two-particle state, accessible 
via quantum state tomography based on current and current correlation measurements, 
we establish that the entanglement of the reduced state constitute a lower bound for 
the entanglement of the projected state. In the two-particle interferometer experiment 
the reduced state is however marginally entangled. Moreover, a finite temperature Bell 
Inequality formulated in terms of currents and current correlators can not be violated in 
the experiment. This shows that an unambiguous demonstration of the entanglement via 
measurements of currents and current correlations requires a reduction of the dephasing 
and the temperature. 
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